Inhomogeneities in the freeze-out of relativistic heavy ion collisions at CERN SPS 
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We study the role of temperature and density inhomogeneities on the freeze-out of relativistic 
heavy ion collisions at CERN SPS. Especially the impact on the particle abundancies is investigated. 
The quality of the fits to the measured particle ratios in 158 AGeV Pb+Pb collisions significantly 
improves as compared to a homogeneous model. 



One of the key motivations for the heavy ion programs at 
GSI, CERN and BNL is to shed light on the QCD phase 
diagram. More specifically, the aim is to gain a deeper 
understanding of the physics of the different phases of 
QCD matter and of the characteristics of the deconfinc- 
ment and chiral phase transition (cf. Fig. 1). 
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FIG. 1: Different phases of strongly interacting matter in 
the T-p plane. Relativistic heavy ion experiments heat and 
compress ordinary nuclei to high temperatures and densities. 
Depending on the bombarding energy, the transition to chi- 
rally restored quark- and gluon matter happens in different 
temperature and density regions, where the characteristics of 
the phase transition are expected to be different. Taken from 
[!]■ 

The current picture of the QCD phase diagram is as 
follows: At vanishing chemical potential (/i = 0) finite 
temperature lattice QCD calculations find a rapid but 
smooth crossover (see e.g. [2]). At large /i one has to 
rely on model calculations, since lattice QCD calcula- 
tions encounter the fermion sign problem. However, sev- 
eral different model calculations (see [3] for a summary) 
suggest a first order phase transition. Combining these 
two results, the line of first order phase transitions origi- 
nating at the T = axis cannot end at /1 = but at some 
point in the (T c , ji c ) plane with finite /x. At this endpoint 
a second order phase transition is expected. For chem- 
ical potentials smaller than /i c a crossover occurs. This 
picture is also supported by different extrapolations of 
lattice QCD to finite chemical potential [4, 5]. 

E.g., in [5] it was found that a line of first-order phase 



transitions in the (/^s,T) plane ends in a critical point 
at T fa 160 MeV, \x B ps 360 MeV (cf. Fig. 2). 
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FIG. 2: Quark-hadron phase transition as obtained from lat- 
tice calculations (at finite chemical potential). Taken from 
[5]. 

Heavy-ion collisions at high energies are hoped to be 
able to detect that critical point and verify this picture. 
For high enough energies the phase transition/crossover 
line should be crossed - the critical energy density is ex- 
pected to be reached at intermediate SPS energies or the 
new GSI facility. Since the system takes different paths 
in the T — fx or T — p plane for different bombarding en- 
ergies, it is hoped that by varying the beam energy, one 
can "switch" between the regimes of first-order transition 
and cross over, respectively (cf. Fig. 1). 

But how do we know whether the system passed 
through a first order phase transition, a crossover, or if 
no phase change at all occured? 

Using a non-equilibrium hydrodynamical simulation it 
was shown in [6] that the expanding fluid develops sig- 
nificant inhomogeneities, if a first order phase transition 
is crossed. These inhomogeneities should also be present 
on the decoupling surface of the hadrons. 

Fig. 3 shows the trajectory of the system within the 
phase diagram for different initial conditions. Depend- 
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FIG. 3: Evolution of the average fluid energy through a 
crossover (dots), and a weak (dashes) and strong (dash-dots) 
first order phase transition. The fat dots indicate time inter- 
vals of « 1.5 fm/c. Taken from [6]. 



ing on these initial conditions, the system either evolves 
smoothly through a crossover or enters the region corre- 
sponding to phase coexistence in the equilibrium phase 
diagram and thus undergoes a first order phase tran- 
sition. The resulting RMS fluctuations of the baryon 
density are shown in figure 4. It can be seen that the 
amplitude of the density contrast is substantially larger 
for a strong first order transition (initial energy density 



1.4eo) than for a crossover (e cq = 2.9eo 
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FIG. 4: RMS fluctuation of the baryon density with initial 
conditions chosen such that the system undergoes a crossover 
(dots), weak (dashes) and strong (dash-dots) first order tran- 
sition as a function of the average baryon density. The fat 
dots indicate time intervals of m 1.5 fm/c. Taken from [6]. 

That means, that if the particles decouple shortly 
after the expansion trajectory crosses the line of first 
order transitions one may expect a rather inhomoge- 
neous (energy-) density distribution on the freeze-out sur- 
face (similar, say, to the CMB photon decoupling surface 
observed by WMAP [7], see Fig. 5). 




FIG. 5: CMB photon decoupling surface observed by WMAP. 
Dark regions correspond to cooler and brighter regions 
corrspond to warmer spots. Taken from [7]. 



On the other hand, if the system expands through a 
crossover transition, as expected for collisions at very 
high energies (/is ~ 0, cf. Fig. 1) it may cool smoothly 
from high to low T and so pressure gradients tend to wash 
out density inhomogeneities. Similarly, in the absence of 
phase-transition induced non-equilibrium effects, the pre- 
dicted initial-state density inhomogeneities [8] should be 
strongly damped. 

Unfortunately, if the scale of the inhomogeneities is 
much smaller than the decoupling volume then they can 
not be resolved individually, nor will they give rise to 
large event-by-event fluctuations. However, because of 
the nonlinear dependence of the hadron densities on T 
and hb they should nevertheless reflect in the event- 
averaged abundances. 

Particle production in relativistic heavy ion collisions 
has been investigated in several works (see for example 
[9]) using (homogeneous) thermodynamical equilibrium 
calculations. In addition, e.g. in [10, 13] extensions 
accounting for strange- and light quark non-equlibrium 
were considered and for example in [11] the role of in- 
medium masses was discussed. Here, we attempt to check 
whether the experimental data show any signs of inho- 
mogeneities on the freeze-out surface. To do this we in- 
vestigate an inhomogeneous fireball at (chemical) decou- 
pling. Perhaps the simplest possible ansatz is to employ 
the grand canonical ensemble and - in extension to the 
homegeneous models - assume that the intensive vari- 
ables T and \jlb are distributed according to a Gaussian 
[12]. This avoids reference to any particular dynamical 
model for the formation and the distribution of density 
perturbations on the freeze-out surface. Also, in this sim- 
ple model we do not need to specify the probability dis- 
tribution of volumes V. Then, the average density of 
species i is computed as 

Pi (T,~p B ,5T,8(j, B ) = (1) 

oo 

dT P{T;T,5T) \ d[i B P^b^b^^b) pi(T, fX B ) , 



with pi(T,fiB) the actual "local" density of species i, 
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the distribution 



and with P(x;x, Sx) ~ exp ^— ^ ^ 

of temperatures and chemical potentials on the freeze- 
out surface. Feeding from (strong or weak) decays is 
included by replacing J> i — > Jj i + B^ . The implicit sum 
over j ^ i runs over all unstable hadron species, with Bij 
the branching ratio for the decay j — > i. For the present 
analysis we computed the densities pi(T, ps) in the ideal- 
gas approximation. The resonances are included up to 
1.5 GeV in mass for the mesons and up to 2 GeV for the 
baryons. The finite widths of the resonances were not 
taken into account and unknown branching ratios in the 
particle data book were excluded from the feeding. Fur- 
thermore, we use a four-dimensional table with 5 MeV 
steps in T and ST and 10 MeV steps in p, and Sp. This 
finite grid-size of course limits our accuracy in determin- 
ing the best fits. However, our approach should be well 
suited to investigate the qualitative behavior of the pa- 
rameters ST and Sps and whether they can significantly 
improve the agreement with the experimental data. 

The data used in our analysis are the particle multi- 
plicities measured by the NA49 collaboration in ^snn = 
17.3 GeV Pb+Pb collisions at CERN SPS. We use midra- 
pidity and in data, both as compiled in [13]. 

Using these data, we perform a \ 2 fit- I.e., we deter- 
mine the minimal value of 
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where r^ xp , r™ odel denote the experimentally measured 
and the calculated particle ratios, respectively, and of 
is the experimental error. We compare two cases: On 
the on hand the homogeneous fit, where the values of 
ST and Sp are set to zero, and on the other hand the 
inhomogeneous fit, where we allow for finite values of ST 
and 5 p. 

We find that the fits improve (lower x 2 /dof) substan- 
tially if ST, Sps are not forced to zero. Table I shows 
the resulting best fits, with and without finite widths of 
the T and \xb distributions. 





T p B ST Sp B X 2 / dof 


SPS-158 
(mid) 


155 ± 5 200 ± 10 40.4/8 
105 ± 5 230 ± 15 35 ± 5 80 ± 40 11.2/6 


SPS-158 

(477) 


145 ±5 210 ± 15 40.0/11 
100 ± 5 260 ± 15 30 ± 5 190 ± 35 5.7/9 



TABLE I: Best fit parameters and xV do / to SPS 158 AGeV 
midrapidity and An data, measured by the NA49 collabora- 
tion. The cases ST — Sp — corresponds to the homogeneous 
freeze-out model and those with finite values for ST, Sp cor- 
respond to the inhomogeneous model. 

As can be seen, for 47r data the resulting best fit 
X 2 /dof values are approximately 3.6 for the homoge- 
neous fit and 0.63 for the inhomogeneous fit. For midra- 
pidity data we obtain \ 2 /dof « 5.1 for the homogeneous 



and \ 2 /dof w 1.9 for the inhomogeneous case. I.e., for 
both data sets the \ 2 P er degree of freedom is consid- 
erably reduced by allowing for inhomogeneities or finite 
widths of the T- and /i-distributions, respectively. In 
other words: for midrapidity as well as for An data ST and 
Sp represent significant paramters. Error estimates for 
the parameters (confidence intervals) are obtained from 
the projection of the regions in parameter space defined 
by X 2 < Xmin + 1 on t° each axis. This corresponds to a 
confidence level of 68.3% if the errors are normally dis- 
tributed. As shown in table I, within this error estimate, 
the best-fit values for ST and Sp are significantly greater 
than 0. 

The resulting particle ratios for the different fits are 
compared to the experimental An data in figures 6 and 7. 
As can be seen especially in figure 7, significant improve- 
ment compared to the homogeneous freeze-out model is 
obtained for the Kaons, the mulit-strange baryons, but 
also for the 0, which couples only to temperature fluctu- 
ations. 
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FIG. 6: Particle ratios (An, compiled in [13]) as measured by 
the NA49 collaboration in SPS 158 AGeV Pb+Pb collisions 
compared to the homogeneous fit (ST = Sfj, = 0) and the 
inhomogeneous fit (ST, S/j, free parameters). 



1.6 
1.5 
1.4 
■ 1.3 
c£ 1.2 
~^ 1.1 
1.0 
=5 0.9 
0.8 
0.7 
0.6 



A 



9 data 

I I inhom fit y /dof =0.63 
A horn fit, x 2 / dof = 3 - 64 




□ 



vr" K + K" K° 



A A S" s o n 



FIG. 7: Fitted particle densities divided by the corresponding 
measured 4n particle density (NA49 collaboration, Pb+Pb @ 
SPS 158 AGeV, compiled in [13]) for the homogeneous (ST = 
S/j, — 0) and the inhomogeneous fit (ST, Sfi free parameters). 
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As can be seen from table I, the inhomogeneous fits re- 
turn significantly lower mean temperature T. However, 
these do not correspond to the "mean" emission tem- 
perature of the particles. The actual particle emission 
distribution is obtained by folding the assumed Gaussian 
(T, /i)-freeze-out distribution (dashed line in Fig. 8 and 
9) with the ideal gas density distribution for a given par- 
ticle species. The resulting normalized probability distri- 
butions read: 

D z (T,T,JI B ,5T,6 f iB)= (3) 

OO 

P(T;T,ST) J rf^i B P( / u B ;7Z B ,(5^ B ) pi{T,pb B ) 

— OO 

Pi (T,~p B ,5T 7 6n B ) 
Di(fj, B ,T,JI B ,ST,SfXB) = (4) 

OO 

P(Vb;~Pb^Vb) J dfi B P(T;T,5T) pi(T,fi B ) 

— OO 

Pi (T,~P b ,$ t ,^b) 
They are shown in figure 8 and 9, respectively. As 
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FIG. 8: Relative particle density versus temperature for 
ir,K,p and p which are obtained from folding the underly- 
ing Gaussian (T, /i)-freeze-out distribution (dashed line) with 
the ideal-gas particle density distribution. 

expected from the temperature distribution we observe 
that the particle emission distributions are shifted to- 
wards higher temperatures. How much the distribution 
is shifted depends on the mass and the chemical poten- 
tial of the corresponding particle species. The probabil- 
ity distributions of the different particle species versus 
chemical potential are shifted to larger chemical poten- 
tials for the particles and to lower chemical potentials for 
the antiparticles. 

Thus, from the finite widths of the Gaussian, differ- 
ent particle emission distributions, with different peaks 
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FIG. 9: Relative particle density versus chemical potential for 
7r, K,p and p. 



for different particle species, result. The corresponding 
means of the distributions can be evaluated as 



< T >i 



< Vb >i 



dT T Di(T, T,~p B ,5T, 8p B ), 



(5) 



/ ' dp B Ms Di{p B ,T,p B ,5T,5p B )- (6) 

The resulting means of the distributions for our anal- 
ysis of the SPS 158 AGeV data are shown in table 11. 



SPS 158 


p p K+ K~ n Q 


< T > [MeV] (mid) 

< hb > [MeV] (mid) 


157 170 152 150 164 180 
268 191 237 222 234 225 


< T > [MeV] (4tt) 

< p-B > [MeV] (4tt) 


136 153 140 139 151 165 
487 22 306 213 277 206 



TABLE II: Mean temperature and chemical potential of var- 
ious particle species for the inhomogeneous freeze-out. 

Again it can be seen that the antiparticles are mainly 
emitted from regions with small chemical potentials while 
the particles mainly originate from high chemical poten- 
tial regions. The mean emission temperatures are in the 
range of 150-180 MeV for the midrapidity data and of 
136-165 for the 47r data. One sees that the antiparticles 
emerge in general from hotter regions than the particles. 

In conclusion, we find that allowing for inhomo- 
geneities in the freeze-out temperature and chemical po- 
tential in an ideal gas description of particle produc- 
tion in heavy ion collisions, significantly improves (lower 
X 2 /dof) the description of experimental data at SPS 158 
AGeV. ft follows that the bulk of the particles orig- 
inates from different density and temperature regions 
than the corresponding anti-particles. Hence, our results 
suggest that the decoupling surface might not be very 
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well "stirred" . Furthermore, inhomogencities appear to 
cure some dcficcncics of homogeneous freeze-out models 
and they might represent a potential variable to connect 
the measured particle abundances to the course of the 
expansion of the system. The investigation of lower SPS 
and RHIC energies within our model is under way [14]. 
Furthermore, in a future comprehensive analysis, the in- 
homogeneities should be generated within a dynamical 
description. 
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